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1 Introduction 

For various reasons as discussed for instance in [1], the selection of actuator and sensor 
positions is still ad hoc. This is especially true for flexible structures where many candidate 
configurations can exist. This study is an attempt to make the selection process more 
methodical. 

One approach to actuator and sensor placement is to optimize a closed loop performance 
metric directly by selecting the actuators, sensors, and controller gains simultaneously. This 
direct approach makes sense if the desired closed loop performance is well defined. Since 
the individual actuator and sensor contributions to the closed loop performance metric is 
complex, the solution strategy usually employs non linear programming with many design 
and numerical iterations. A second approach is to select actuators and/or sensors based 
on open loop properties so that closed loop performance is indirectly optimized. Since the 
individual sensor and actuator contributions to the open loop metric is simple, nonlinear 
optimization is usually not needed. This approach will suggest efficient actuator and sensor 
configurations for any type of control law. The method suggested in this study falls into the 
latter class of approaches. 

In this note, the use of Hankel singular value (HSV) formula as an actuator placement 
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metric [2] is extended to flexible structures in discrete time [3]. A main novelty introduced 
is that the ambiguity in the weighting of the principal modes is addressed by incorporating 
the general disturbance rejection goal into the actuator and sensor placement formulation. 
Optimal actuator and sensor placement is considered for the purpose of designing control 
laws for the general disturbance rejection problem. This apparent restriction to disturbance 
rejection problem is not too restrictive sine# it is well known from modern multivariable 
control theory that stability and even robustness requirements can be transformed to this 
form with appropriate weighting. Simulation results demonstrates that the improvement 
in closed loop performance is independent of the type of controller used since open loop 
properties have been improved. 


2 Actuators and Sensors for Disturbance Rejection 

Figure f show a schematic of the disturbance rejection problem where the inputs to the plant 



Figure f: General disturbance rejection problem 


consist of two vectors and similarly the outputs consist of two vectors. As a necessary starting 
point in a disturbance rejection problem, the disturbance source (or inputs) and the output 
response defining the performance to be optimized is usually defined by the requirements of 
the control problem. 

The disturbance rejection viewpoint taken in this study is similar to the mode selection 
viewpoint proposed earlier in the context of model reduction for flexible structures [4]. The 
above formulation is also consistent with the general multivariable control design [5] where 
the plant P = [P n , P r2 ; P 2 i, P 22 ] is assumed given and the problem is to design a stabilizing 
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and realizable controller, K to minimize a suitable norm of the closed loop transfer function 
matrix, Fi(P } K). The basic difficulty in the actuator /sensor placement problem is that only 
Pu is given. 

In this note, HSV are used to construct a metric that quantifies the degree of control- 
lability and observability for a given set of sensor and actuator configuration. Although 
the use of HSV to analyze the degree of controllability and observability of a linear system 
is well established, especially in model reduction applications [4, 6, 7], the approximate 
decomposition of the HSV with multiple sensors and actuators in terms of the HSV of all 
combinations of sensor and actuator pair is new. This result significantly simplifies the design 
problem of selecting the most effective set of sensors and actuators for flexible structures. 
The main novelty of the placement strategy introduced in this note is that the HSV from 
the disturbance to the performance outputs are used to weigh the HSV between candidate 
actuators to sensor sets. This is possible because in both cases, individual HSV directly 
corresponds to individual structural modes which is unique to flexible structures. 


3 Decomposition of HSV 


Given the quadruple [A z , B Z ,C Z , D z ) of a discrete LTI state space matrix of a flexible 
structure, assumed to be lightly damped with distinct eigenvalues, let (A, B } C\ D) denote 
the 2 by 2 block diagonal form whose zth block is 


Re(^i) -Im(^) 
Im (z t ) Re(zi) 


( 1 ) 


Let (zi 7 Vi) denote the zth eigenvalue and eigenvector pair of A z and T denote the sampling 
period. The steady-state discrete time controllability grammian, W Coo , and observability 
grammian, W 0oo , satisfies Sylvester equations. The triple (A } B } C) is internally balanced 
if its grammians are equal and diagonal [7], i.e. , W Coo = W 0oo = T 2 . The nonnegative real 
diagonal elements of T are called the Hankel singular values of the system. 

Due to the diagonal dominance property of the discrete controllability and observability 
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grammian for flexible structures, the square of the zth HSV is given [3] by 


trace BB T ..trace C T C 

,4 rsj L J_U L J U 

8 “ (4 s t ry 


( 2 ) 


where Si = — ^Re^nzi) and the subscripts ii denotes the zth 2 by 2 block of matrices formed 
from the inputs and outputs. It is shown [3] that the approximate formula is quite accurate 
up to frequencies near 90 % of Nyquist frequency. 

For p actuators and q sensors, the input and output matrices consist of p columns and 
q rows, respectively 


B z = [B Z1 ,..., B Zp ], 


cj = [c%. 


c 1 

•5 ^ Z 


(3) 


Partition the state transformation matrix as follows 


V=[r 1 ,...,r„]; = [if , . . . , if 


(4) 


where rq = [Re(ry) — Im(ry)], and if are 2 n by 2 matrices made up of components of zth 
right and left eigenvectors. The approximate square of the zth HSV in Eq.(2) can be written 
in terms of individual actuator and sensor 


4 „ e u ffj • n=i gi 

li = - 


(4 S t Tf 


j = 1 k = 1 


where 


ffk ~ \\kB Zk \\l; g 2 lk — 1 1 Zk 1 1 2 


(5) 


( 6 ) 


7 8 4 (j, k ) denotes the squared zth HSV for the jth input and kt h output. 


The importance of Eq.(5) lies in the decomposition of the HSV of the multivariable 
flexible structure in terms of the sum of approximate HSV of each actuator and sensor 
pair. The contribution of each pair appears in a very convenient form from their placement 
standpoint. The exact HSV are not decomposable in general, i.e., 


it = ME E ^ EE 

j = 1 k = 1 j = 1 k = 1 


(7) 
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4 Placement Metric 


Let denote the HSV from the disturbance to the performance outputs (see Figure f) 

r WZ diag iwzn ) (8) 

Also let Y uy and Y uy denote the HSV for p actuators to q sensors and the HSV from a 
reference set of actuators to a reference set of sensors respectively 


^uy = diag (Aim/i ■>•••! luy n ) 
r uy = diag (7 uyi , . . . , 7 -ujm ) 

Dehne the actuator and sensor placement metric as a weighted sum 


71 /y 4 

7 Ay ’uyi y 

° / -j -A 


8 = 1 h-UJ/, 


4 iwzi 


q p 

= E E 4* 

j=l k = 1 


( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 


where contribution from jth sensor and kt h actuator pair is 

7=1 \ luy t / 

It is necessary to introduce the above weights because the basic idea is to improve the joint 
controllability and observability for multiple modes and not all modes are equally important 
physically. HSV can vary by orders of magnitude and therefore a few physically irrelevent 
modes could dominate an unweighted scalar metric. Hence, the normalizing weight, 7^., is 
introduced to make each HSV equally important assuming no a priori physical knowledge 
of individual modes. An example of a normalizing weight is when all actuators and sensors 
are included. After the normalization, the relative importance of each mode in the control 
problem can be incorporated by using the weight, 7^ . Modes that are important in the 
disturbance rejection performance will give larger weights. Note that Jjk is basically a 
weighted trace of T uy (j, k) squared. 

For actuator (sensor) placement only, the set of q sensors (p actuators) are assumed 
fixed. The above metric simplifies to 


J = E -C E f 


( 13 ) 


k = 1 


V'=i 


5 



where the contribution of the kt h actuator (jth sensor) is 


J a k ct = £ /, 


ik w i 


j; en = 


2 2 


8 = 1 


8 = 1 


The weighting factor for the zth mode is 

l ^ i 


wj = 


A 

I WZ, 


(i6,T) 2 li„ JUS 


,(4<W) 2 7 4 


-A P 

/81^ r 2 

/ . 1 8/8 


«8/t fc = l 


(14) 


(15) 


Note that denotes the influence of the kt h actuator (jth sensor) on the zth mode. 

Additional factors that contribute to the above weighting factor include damping, significance 
with respect to disturbance rejection, and the reference HSV, for the zth mode. 


5 Placement Strategy 

Table 1 show the individual contributions for all candidate sensors and actuators. The sum 
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Table 1: Actuator and sensor selection table 


of each row is given by the last column while the sum of each column is given by the last row. 
The total metric for all sensors and actuators are given by the sum of the last row or column. 
The costs J/ ct ( Jj £n ) represents the contribution of the kth actuator (jth sensor) using all 
given sensors (actuators). The goal is to maximize the metric J using the least number of 
actuators and/or sensors. Note that the advantage of the metric applied to flexible structure 
is that the contribution of each actuator and sensor appears linearly and independently. 

For simultaneous actuator and sensor selection problem, we need to maximize the linear 
sum of non-negative numbers 

J = (W) 

j = 1 k= 1 
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Clearly, a table of nonnegative numbers need to be examined such that J can be approached 
using the fewest number of actuators and sensors. For actuator (or sensor) only selection 
problem, a bar chart of 

j = = mi 

k = 1 j = l 

will indicate the best selection. Actuators (sensors) with small values of J^ ct (J s - en ) can be 
removed as the least significant ones. 

6 Concluding Remarks 

The metric and methodology introduced in this paper appears promising in its potential 
development as a generic tool for actuator and sensor selection for feedback control design 
of flexible structures. This optimism is based on the advantages evident in using the 
approximate HSV formula for actuator /sensor selection: HSV of individual pairs can be 
summed; and it is analytically and computationally simple and avoids direct numerical 
optimization; and perhaps most importantly, it gives physical insight. 

Preliminary simulation results based on i7 2 and H control law designs indicate that 
closed loop performance can be improved significantly and the performance improvement is 
independent of the particular type of control law. This consistent improvement is clearly the 
result of judiciously improved controllability and observability in the open loop system. How- 
ever, the particular form specified for the disturbance rejection performance requirements 
strongly affect optimal actuator selection. 
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